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1. INTRODUCTION 

The majority of electric power distribution lines in India are overhead lines. These lines are subject 
to disturbance such as an energized conductor falls to the ground or connected to high impedance objects due 
to dissimilar environments. This condition is called High impedance fault (HIF), and it draws little fault 
current which is less compared to the threshold values of conventional protection means. Hence, some of 
these HIFs may not be detected. Further, threshold values cannot set as lower values, since the normal load 
distribution will lead to nuisance tripping. Besides, HIFs are accompanied by electric arcs. These arcs 
produce an asymmetric, unpredictable, and random current signal. Most of the electric power distribution 
networks (EPDN) are in close proximity to the thickly populated area. If EPDN is operated with HIF leads to 
endangering the human lives and their properties or material goods. Therefore, with the goal of increasing the 
performance and the safety of EPDN, several papers have been published in the past few decades to solve 
HIF problem is well documented in [1]. The redefining process of detection yet to be completed due to an 
increase in non-linear loads (NLLs). 

A method based on harmonic content presented in [2, 3]. The setback of such method is to set 
threshold values which reduces the performance of the method. Time-frequency analysis [4, 5] produces 
better performance in HIF detection. The modified Fast Fourier Transform (FFT) method employs the 
relative relation between the third, fifth and seventh harmonic current in [6]. The presence of NLLs taken 
into account and reliable detection shown in the system studied. However, the percentage of false tripping 
and computational work is an obstacle for practical implementations. Time-domain based method called 
mathematical morphology [7, 8] is used for detection of HIF, which is good for a balanced system. When the 
system is unbalanced, these methods show less performance. 

Wavelet Transform (WT) has been used in HIF detection. WT is analyzing the signal with 
frequency component and their position in time. More than 20 years such methods used in protection 
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applications. The algorithm in [9] uses dynamic features extracted by stationary WT and a support vector 
machine based decision-making system. Arcing currents related to different high impedance surface are 
generated in the laboratory setup, and these signals are decomposed using Discrete Wavelet Transform 
(DWT) [10]. The authors used DWT to examine high- and low-frequency voltage components at different 
points [11] and an average of absolute difference of extracted voltage signal [12] of EPDN. In [13] voltage 
and current signals are used to detect and locate HIF respectively. This method uses the absolute sum of 
high-frequency components for both detection and location. Two different HIF detection scheme are 
discussed in [14], both utilizing WT decomposition. In the first method, principal component analysis (PCA) 
for feature reduction and Neural Network (NN) for feature classification. In the second method, a Genetic 
Algorithm (GA) for feature reduction and Bayes clasifier for classification. The change in current signal 
created by HIF and other transient events has been used in [15] with NLLs in the system. A DWT used to 
decompose the current signals and extracts features to train NNs. An evolving NN [16] and continually 
online trained NN method [17] was shown as the proper approach to HIF detection since it is a time-varying 
problem. A DWT with data mining based classification in [18]. 

These research articles disclose useful properties and different detection methods of HIF. However, 
the majority of the above-mentioned articles fail to consider NLLs except in [6, 15, 18]. These three articles 
are showing less variation in NLLs. Since NLLs such as television, computer, fluorescent lamps, etc., 
constantly increasing year by year in the distribution grids and NLLs and HIFs characteristics have closely 
resembled each other, which will make the existing methods less effective. Hence, an enhanced detection 
method for HIFs with a massive variation of NLLs is Proposed. 

The rest of the paper has been prepared as follows. Section 2 explains the test system and also the 
characteristics of HIF model engaged for simulation. The proposed methodology has been thoroughly 
discussed in Section 3, including the decomposition of a signal, feature extraction, feature selection, and 
classification. Section 4 has results and discussion and the paper concluded with the main highlights of the 
work in section 5. 


2. TEST SYSTEM 
2.1. System studied 

The proposed method is verified through an actual EPDN in India (Chennai) shown in Figure 1 was 
modeled with a sim-power-system block set (MATLAB). The system data are given in [18]. Distribution 
lines are modeled as lumped parameters. Busbar input is modeled by Thevenin's voltage and 
equivalent impedance. 


2.2. HIF model and characteristics 

The HIF Emanuel’s model [4] shown in Figure 2 used in this work. The HIF current signal during 
different parameters of the model is shown in Figure 3, which shows the non-linear behaviour, asymmetry in 
the signal, fewer current and random behavior of fault current. This is due to low- and high-frequency 
content of the signal is changing with time and hence the signal is termed as a non-stationary signal. 


2.3. Test conditions 
For the classification HIF and non-HIF (NHIF), a total of 540 HIF and 460 NHIF cases are 
described in Table 1. 


3. PROPOSED METHODOLOGY 

The proposed HIF detection approach is shown in Figure 4. The process starts with the simulation of 
EPDN and accessing current signal. The DWT decomposes the current signal and features such as energy, 
entropy, skewness, sum, standard deviation, and kurtosis are extracted. These features are training two layer 
feed forward NNs to classify the case belongs to HIF or NHIF. 


3.1. Signal decomposition using wavelet transform 

Wavelet is an effective tool to analyze non-stationary signal with the capability of multiple 
resolutions in time and frequency. DWT is a flexible signal processing method broadly used in power system 
engineering. DWT has been successful in stability analysis of power system engineering. The mathematical 
equation for DWT represented in (1). 


n-k2™ 


DWT mn) = RUX(K) WA) (1) 
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where x(k) is the discrete signal in terms of coefficients. 
y(.) is mother wavelet, m and n are time scale parameters. 
k is discrete time of coefficients. 

2™ is scaling parameter. 

k 2™ is shifting parameter. 


1 . ag 
— is the energy normalization component. 
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Figure 1. EPDN studied Figure 2. HIF model 
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Figure 3. HIF current waveform 
(a) RI=10000 Q, R2=9000 O,V1=5000 V,V2= 5050 V. 
(b) R1= 150 Q, R2= 160 Q, V1=2000 V, V2=2050 V. 


Low pass filters (LPF) and high pass filters (HPF) is used in DWT to decompose a signal into a low- 
frequency and high-frequency component. The decomposition process starts by passing a signal through LPF 
and HPF. The LPF produces an approximate coefficient which has high scaled and low-frequency 
decomposition. In contrast, the HPF produces detailed coefficient which has low scaled and high-frequency 
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decomposition. At once the first level of decomposition is completed, the sampling frequency is reduced by 
half of its value. To calculate the next DWT level, LPF output (approximation) is decomposed. The 
Daubechies basis db4, five level decomposition, and the sampling rate of 25.6 kHz with the data window 
length of three cycles are chosen in this work. 


Table 1. Test conditions 


Event Simulated conditions Total 

HIF Resistances varied from 150 Q to 12 kQ and DC voltage from 1.5 kV to 10.5 kV 540 
randomly 

Load switching Changes in load level: 0% - 25%, 25% - 50%, 50%-75%,75%-100% and 100% -110% 90 
both forward and reverse conditions. 

NLL switching The ratio of linear to non-linear load changes 1 to 0.5 in steps of 0.05 both in forward and 216 
reverse conditions. 

Capacitor Switching conditions on and off with load variation of 25%, 50%, 75 %, 100% and 125%. 90 

Switching 

Inrush current No load switching of 33 /11 kV transformer 64 
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Figure 4. The proposed HIF detection approach 


The decomposition of phase a current due to NLLs and capacitance switching at inception time of 
0.065 s is shown in Figure 5 and Figure 6. It is found that the behavior of both NLL and capacitance has low 
and high-frequency parts in all wavelet detail coefficient DS to D1. The decomposition of phase a current 
during HIF is shown in Figure 7. A small magnitude of change in frequency component is present since the 
initiation of HIF in detail coefficient D5 alone. 


3.2. Selection of wavelet coefficient 

From the DWT analysis of the current signal, there were one approximation coefficient (A5) and 
five detailed coefficients (dl-d5) acquired. The comparison of different wavelet detail and approximation 
coefficient has been made to select wavelet coefficient for HIF Detection. It is observed that detail 
coefficients d4 and d5 offers diverse characteristics during different power system events and hence selected 
for feature extraction. 


3.3. Features using wavelet coefficients d4 and d5 
Six different features were estimated using wavelet coefficients d4 and d5 are described below: 


a. Standard deviation: It represents the deviation of a signal from its mean. 
b. | Energy: The total energy content of the current signal. 


Energy = JN i? (t) dt (2) 


c. Kurtosis: Bigger kurtosis point represents more outlier in the signal. 





n(n+1) (x—x)* 3(n-1)? 
Po aao 6) 


Kurtosis = ESS (n-2)(n-3) 


o* 
o- Standard deviation, x — Mean, n — no of sample data 


d. Skewness: It is a measure of the irregularity of the probability distribution about its mean. 
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Figure 5. DWT of current signal under non-linear load switching 
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Figure 6. DWT of current signal under capacitance switching 


Q 
gz 
5 0.2 
= 0 
< 
-0.2 
0 0.05 0.1 0.15 0.2 
Detail D2 
0.05 
0 
-0.05 
-0.1 
0 0.05 0.1 0.15 0.2 
Time (s) 
j Approximation A5 
0 
ait 
0 0.05 0.1 0.15 0.2 
Detail D4 
0.04 
E 0.02 
: 
G 0 
a -0.02 
-0.04 
0.05 0.1 0.15 0.2 
Detail D2 
0.02 
0 
-0.02 
0.05 0.1 0.15 0.2 
Time (s) 


Detail D5 


0.05 0.1 
Detail D3 


Amplitude 
= 


0.05 0.15 0.2 


0.1 
Detail D1 


0.05 0.1 0.15 0.2 
Time (s) 


Figure 7. DWT of current signal under HIF 
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e. Entropy: It gives complexity of the signal. 

Entropy = — } pj log pj (5) 

Where pjis an energy probability distribution of WT detail coefficient (d4 and d5) 
f. Sum: Sum of all points of signal. 


3.4. NN implementation 

The NN is a powerful tool used in pattern recognition and classification [19]. Consequently, the 
proposed work uses NN to learn the input and output relationship from feature input vector. As mentioned in 
the feature extraction process, the input to NN comprises 12 elements (6 features of detail coefficients d4 
and d5). The output layer has two neurons 1.e. either HIF or NHIF and 10 hidden layers. A two-layer feed- 
forward NN has been used in this work as shown in Figure 8. The sigmoid activation function used for 
hidden layer and softmax for the output layer. The NN is trained with scaled conjugate back-propagation that 
updates weight and bias values. Also, the number of epochs was determined by experimentation. The number 
of epochs for the training was set around 1000 and a mean squared error rate of e™® was used [20]. The 
dataset is randomly divided into three parts: a training data set of 70%, validation of 15% dataset, and the 
remaining 15% for testing. 


Hidden Output 





Figure 8. NN model for HIF detection 


4. RESULTS AND DISCUSSIONS 

A set of 1000 features obtained by changing the simulation conditions as given in Table 1. Thus a 
wide range of studies has been carried out to demonstrate the proposed method. Conversely, it is not possible 
to report all information because of space constraint. The proposed method is evaluated through the 
following three parameters: 

- Dependability: Predicted HIF against total HIF conditions. 

- Security: Predicted Non-HIF against total non-HIF conditions. 

- Accuracy: Actual predicted against the total number of conditions considered. 

- Speed: One cycle power frequency current / No of cycles that it takes to detect the fault. 

The confusion matrix of the NN during training, validation, and testing are presented in Figure 9. 
During training, all 370 HIF test cases are correctly classified to produce dependability of 100%. Two NHIF 
test case is misclassified as HIF to produce 99.4% security with the overall accuracy of 99.7%. During 
validation of 150 cases, dependability observed as 100%, and 97% of security with the overall accuracy of 
98.7%. Finally, during testing of 150 cases, all HIF cases predicted correctly and only one NHIF cases are 
classified as HIF and produce dependability, security, and accuracy as 100%, 98.4%, and 99.3% respectively. 
It is very clear from the confusion matrix that the method has an excellent detection rate. Nevertheless, there 
was a huge amount of non-linear load variation, all the HIF cases are detected. 

Power system more vulnerable to noise, it is necessary to investigate the proposed method with 
different Signal to Noise Ratio (SNR). Table 2 compares the performance under noisy atmosphere. It noted 
that the performance is similar to the one observed in normal case for 10 dB of SNR on the signal. For 20 dB 
of SNR, HIF detection rate alone same. However, the accuracy and security are reduced to 98% and 95.16% 
respectively. Additional SNR reduces all the three performance indices. Though the accuracy of 20 dB is 
less, the HIF detection rate is excellent, and hence the proposed method is suitable for SNR of 20 or less. 
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Figure 9. The confusion matrix of NNs 


Table 3 shows the impact of the sampling rate of the proposed work. There was substantial 
upgrading demonstrated by the proposed methodology with sampling rate from 256 to 512. The detection 
rate diminishes as the sampling rate value decreases. Hence, the proposed scheme is better with 512 
samples per cycle. 

Finally, the performance of the proposed HIF approach is compared with previously published 
articles as presented in Table 4. It is observed that all three performance indices are better than [4, 7, 10, 14]. 
The speed of detection is compared with the references mentioned in Table 4. The speed of detection is low 
compared to [4] and as good as with other references listed in Table 4. However, it is noted that HIFs are low 
current, during which the distribution network are not in stress. As a result, the HIF detection speed is not a 
critical concern. In other hand, the security index of the HIF detection algorithm is important as well as its 


dependability index. 


Table 2. Impact of SNR on the proposed method 


Performance indices 10 dB 20 dB 30 dB 
Accuracy 99.33% 98.00% 93.33% 
Security 98.38% 95.16% 88.70% 


Dependability 100% 100% 96.59% 


Table 3. Impact of sampling rate on the proposed method 


Performance indices 512 420 256 
Accuracy 99.33% 96.00% 86.66% 
Security 98.38% 95.16% 80.64% 
Dependability 100% 98.88% 90.90% 


Table 4. Comparison of the proposed method with previously published work 
Method Accuracy (%) Security (%) Dependability (%) Speed 


Ref [4] 93.6 81.5 100 1 
Ref [7] 97.3 96.3 98.3 0.04 
Ref [14] 96 100 90 0.25 
Ref [10] - 68 72 0.33 
The proposed method 99.33 98.38 100 0.33 
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5. CONCLUSION 

In this paper, a new methodology based on the WT NNs is presented. The proposed method 
efficiently distinguishing HIF form other normal patterns. The proposed method comprehensively studied an 
actual electric power system with real data using MATLAB with a huge variation in power system operating 
conditions. The method has an exceptional detection rate even under wide variation in NLLs. This indicates 
that the proposed scheme is highly consistent and secure detection for HIFs with a dependability 
index of 100%. 
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